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Abstract 

Measurements of the kinetic Sunyaev-Zel'dovich (kSZ) effect from instruments such as the South 
Pole Telescope (SPT) and the Atacama Cosmology Telescope (ACT) will soon put improved con- 
straints on reionization. Popular models assume that UV photons alone are responsible for reion- 
ization of the intergalactic medium. We explore the effects of a significant contribution of X-rays 
to reionization on the kSZ signal. Because X-rays have a large mean free path through the neu- 
tral intergalactic medium, they introduce partial ionization in between the sharp-edged bubbles 
created by UV photons. This smooth ionization component changes the power spectrum of the 
cosmic microwave background (CMB) temperature anisotropies. We quantify this effect by running 
semi-numerical simulations of reionization. We test a number of different models of reionization 
without X-rays that have varying physical parameters, but which are constrained to have similar 
total optical depths to electron scattering. These are then compared to models with varying levels 
of contribution to reionization from X-rays. We find that models with more than a 10% contribu- 
tion from X-rays produce a significantly lower power spectrum of temperature anisotropies than 
all the UV-only models tested. The expected sensitivity of SPT and ACT may be insufficient 
to distinguish between our models, however, a non-detection of the kSZ signal from the epoch of 
reionization could result from the contribution of X-rays. It will be important for future missions 
with improved sensitivity to consider the impact of X-ray sources on reionization. 
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I. INTRODUCTION 



The epoch of reionization (EoR) is an important milestone in the history of our universe 
There has been much theoretical work predicting the details of reionization, however 
many of the predictions have yet to be verified observationally. An important probe of 
this process is the small scale temperature fluctuations in the cosmic microwave background 
(CMB). In particular, secondary anisotropics are imprinted on top of the primordial CMB 
from the Doppler shift associated with CMB photons scattering off free electrons in between 
the observer and the surface of last scattering at z ~ 1100. This is known as the kinetic 
Sunyaev-Zel'dovich (kSZ) effect. Because only ionized regions can contribute to this effect, 
the power spectrum of kSZ fluctuations is sensitive to details of the reionization history 
00. 

The most popular models of reionization assume that the UV photons from early stars 
drive the ionization of the intergalactic medium (IGM). Due to the short mean free path 
(MFP) of UV photons through the neutral IGM, the resulting topology of reionization is a 
two phase IGM consisting of ionized bubbles surrounded by regions of neutral gas Q. 

Another possibility is that X-ray photons make a significant contribution to reionization 



e.g. |8Hl2|). These X-rays could originate from accretion onto massive black holes in 



galactic nuclei, X-ray binaries, or inverse Compton scattering in supernova remnants [13 
Because the MFP of X-rays through the IGM is much longer than UV photons this will 
result in a more uniform background of ionizing photons. Thus, a significant contribution of 
X-rays can qualitatively change reionization by creating a partially ionized IGM in between 
the bubbles produced by UV photons. 

In this paper we use semi-numerical simulations to gauge the imprint of X-ray sources 
during reionization on the kSZ signal. Because of the smooth ionization caused by X-rays, 
the so called "patchy" component of the kSZ power spectrum can be substantially reduced. 
We find that models with more than a 10% contribution from X-rays produce a significantly 
lower power spectrum of temperature anisotropies than all the UV-only models tested. 

The paper is structured as follows. In §2, we describe our method of estimating the kSZ 
signal including a description of our semi-numerical simulations of reionization with and 
without a contribution from X-rays. In §3, we report our results, which include changes 
in the CMB power spectrum due to the contribution of X-ray sources to reionization. In 
§4 we summarize our main conclusions. Throughout, we assume a ACDM cosmology with 



tt A = 0.73, Vt m = 0.27, n b = 0.0456, h = 0.7, and a s = 0.81 14 



II. METHOD 

A. The Kinetic SZ Effect 

The kSZ effect is the contribution to fluctuations in the CMB temperature associated 
with the Doppler shift of CMB photons Thompson scattered off free electrons moving at 
some radial velocity relative to the observer. The kSZ temperature fluctuation in a particular 
direction is given by 

AT f 

— = a T / dr]e~ T{v) an e v\\, (1) 

where <tt is the Thomson cross section, t{j]) is the optical depth from the observer to the 
position corresponding to conformal time r], a is the cosmic scale factor, n e is the number 
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density of free electrons, and v\\ is the component of the velocity along the line of sight. 

We wish to calculate the impact on the kSZ effect from an X-ray contribution to reioniza- 
tion using semi-numerical simulations of the EoR. These simulations will examine varying 
levels of X-ray contribution to reionization. We generate a number of simulation cubes that 
each span a solid angle of one square degree on the sky (corresponding to a comoving length 
of L 150 — 200 Mpc from z = 6 — 20) and have a resolution of 256 3 separate volume ele- 
ments. In each 3D element of the simulation we calculate the matter over-density, velocity, 
and ionization fraction. Cubes are produced at appropriate cosmic times and stacked along 
the line of sight to give continuous coverage from a redshift of z = 6 to z = 20, corresponding 
to the EoR. We then perform the integral in Eq. (Cp) along each line of sight to produce a 
2D map of CMB temperature fluctuations. 

We compute the density and velocity fields of our simulations using the 21cmFAST soft- 
ware package 15[. Non-linearities are included using leading order perturbation theory 



(Zel'dovich approximation). We expect this to be a reasonable approximation, as the power 
spectra of the density field using this approach agree well with that from detailed hydrody- 



namic simulations at the redshifts and scales relevant for our calculations [151 ] . We generate 
our ionization fields with the model described below. Note that we do not produce data 
cubes at redshifts below z = 6. The kSZ fluctuations from these redshifts will be independent 
of reionization history and we incorporate their effect analytically. 



B. Reionization Model 

To simulate the reionization history we use a semi-numerical approach similar to that 
described in [Til] - The basic assumption in our prescription is that there is a fixed number 
of ionizing photons, A^ on , released into the IGM for each baryon which is incorporated into 
galaxies. From extended Press-Schechter the ory it is possible to calculate the rate at which 
matter will collapse into dark matter halos [16J. In this way one can calculate the ionized 
fraction in a region of a given size and density analytically. For the case without X-rays where 
all ionizing photons are assumed to be UV, the evolution of the ionization in a spherical 
region of comoving radius R with over-density S relative to the cosmic mean is given by [l7| 

dQn, s N ion f dF(M ion ) ^ dF(M min ) \ 

Qr,8 T7 + {I-Qr,6) t: ~ (2) 



dt 0.76 V dt y ' dt 

where the 0.76 accounts for the cosmic mass fraction of hydrogen, Qr^ is the mass averaged 
ionized fraction of hydrogen, and F(M) (which depends implicitly of 5 and R) is the fraction 
of mass which has collapsed into dark matter halos above mass M, 

f (M) = erfc| f 51-* ( '7,°> )■ (3) 



R 



where 5 C = 1.686/ D(z) is the threshold for spherical collapse at z ^> 1 extrapolated to 
redshift zero with linear theory, D(z) is the linear growth factor of density fluctuations, 
5(z = 0) is the linear over-density of the region at the current time, and the cr's are the root 
mean squared density fluctuations on scales corresponding to mass M and radius R today. 
The minimum mass of dark matter halos that can host galaxies is denoted by M min in the 



3 



neutral IGM and by M- lon in the ionized and photo-heated IGM. We parameterize this mass 
in terms of the virial temperature which for z ^> 1 is given by 

/ m \ 1/3 n+ z 



The recombination rate is given by 

B = a B Cr^(l+z) 3 (l + S)Q R)S , (5) 



where ae = 2.6 x 10~ 13 cm 3 /s is the case B recombination coefficient for T = 10 4 K [181 ]. 
C = (n 2 e ) I (n e ) 2 is the dumpiness of the ionized gas in the IGM and is the density of 
hydrogen gas today. 

In order to produce the 3D simulation of the ionization fraction we smooth a random 
realization of the linear density field in our simulation in spheres of varying radius and solve 
Eq. 02]) as a function of the size and over-density at each point. At this stage, regions for 
which Qr^s > 1 are taken to be bubbles and the ionization fraction is set to one while other 
regions are taken to be completely neutral. Next, the equation is solved on the scale of each 
spatial pixel and then used to determine the ionization fraction in regions not encompassed 
by ionization bubbles. 

We use a similar approach to calculate the ionization fraction when X-rays make a sig- 



nificant contribution to reionization. This approach has been adapted from [17|. The key 
simplifying assumption we make is that X-rays have a fixed comoving MFP and that X-ray 
photons uniformly ionize regions of that size. The comoving MFP of X-rays of energy Ex 



through a neutral medium is given by [13j, [19 



1 + zV 2 / E x x '' 



x ~ 180 1— ) UkWj Mpc - (6) 

We assume a constant MFP of A = 20 Mpc, but find that the measured angular power 
spectrum of temperature anisotropics is essentially unchanged if we set A = 50 Mpc. UV 
photons have a very short MFP compared to the size of the spatial pixels in our simulation. 
Thus, when considering a region smaller than the MFP of X-rays one must consider the 
X-ray sources from a larger surrounding region with a radius equal to the X-ray MFP, 
A, while the ionization from UV photons will only depend on the local over-density. We 
capture this behavior by solving a system of 6 coupled ordinary differential equations (ODEs) 
which compute the ionization history (including the relative amounts from X-rays versus UV 
photons) on various scales. As in the UV only case, we smooth the density field on all scales 
to identify ionized bubbles, but in addition simultaneously track the ionization history and 
production of X-ray photons from surrounding regions with radius equal to the X-ray MFP 
(A). The system of 6 ODEs (adapted from [17|) is 

dQn _ dQ^uy dQ^x 



dt dt dt 

dQx _ dQx^uv dQ\,x 
dt dt dt 



(8) 



dQ R ,uv , > N ion ( dF(M ioQ ) dF(M miLII 

K l - A frac J— — Q R — h (1 - Qr) "J: - #3 (9 J 



dt y 11<iW 0.76 V dt v ^ ' dt 
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dt y lldW 0.76 V dt K ^ ' dt 
^^^^ + (l- Qji )!?^)- Bs if «>A (11) 



dt o.76Vi + <wv dt v " 7 dt 



(12) 



dQ A ,x y iV ion / dF(M ion ) , dF(M mm , . „ 

-A frac — — g A — - — + (i-g A ) — - — -±? 6 , (13) 



dt 0.76 V dt v 7 dt 

where Q R and Qa are the mass averaged ionization fractions on scales of R and A respectively. 
Qr,uv is the ionization on scale R due to UV photons and Qr ; x, Qx,x, and Q\,uv are the 
various ionization fractions from either X-rays or UV photons ionization on the two scales 
of interest. The fraction of ionizing photons which are X-rays is denoted by Xf rac . Note that 
the time derivative of the collapsed fraction in each equation is evaluated on the same scale 
as the Q value it is multiplied with. 

The various Bi terms are recombination rates. In the UV case, Eq. (jSJ) was derived 
by assuming that all regions are either completely ionized or neutral. With a significant 
contribution from X-rays this will not be the case. We assume that X-rays uniformly ionize 
the IGM on scale A and UV photons completely ionize bubbles in the IGM around their 
sources. In order to calculate the recombination terms we calculate the volume of the IGM 
which is uniformly ionized, Kinifm-m , and the fraction in ionized bubbles, Vbub- The uniform 
X-ray ionization also allows bubbles to grow larger than they would from the UV photons 
alone; we denote this additional "fringe" volume with Vfringe- These add up to the total 
volume of the spherical regions we consider, V, and are described by 

V huh /V = Q uv (14) 

VMnge/V = ~ (15) 

1 — 

^uniform V ^bub ^fringe; (l 6) 

where either the UV or X-ray ionization must be evaluated on the relevant length scale. It 
follows that various recombination rate terms are given by 

B 3 = a B Cn° H (l + + 5 R )^, (17) 



B B = « u <VJ(l + c)»(H QLv^T^Y (18) 



,0 n i „\3ft i x \ / ^fringe , ^2 ^uniform 

— + y R ,x— —j 

Bi and Bq are the same as B% and B§ except that the appropriate values for 5\ and the vol- 
umes corresponding to scale A must be used. Note that this treatment of the recombination 
rates was not included in previous work. We do not include heating due to X-rays which 
could further suppress the formation of low mass galaxies. This would be straightforward 



to incorporate in the future [17 
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TABLE I: Parameters for the different reionization models used. For the minimum masses of dark 
matter halos which host galaxies, the corresponding virial temperature (assumed to be constant) 
is listed instead of the mass (see Eq. @ in the text). 



Model 


N ion 


C 


T v ; r (M m i n ) 


X v i r (M; on ) 


-^frac 


A 


T 


Fiducial 


60 


2 


10 4 K 


10 5 K 







0.088 


Clump 


355 


20 


10 4 K 


10 5 K 







0.088 


M min 


115 


2 


5 x 10 4 K 


10 5 K 







0.088 


10% X-ray 


53 


2 


10 4 K 


10 5 K 


0.1 


20 Mpc 


0.088 


30% X-ray 


45 


2 


10 4 K 


10 5 K 


0.3 


20 Mpc 


0.088 


50% X-ray 


40 


2 


10 4 K 


10 5 K 


0.5 


20 Mpc 


0.088 



We smooth our density field on all scales and if Qr > 1 the region R is set to be completely 
ionized. After the entire cube has been searched on all relevant scales we solve Eqs. (!71 fT3|) on 
the scale of each spatial pixel to determine the ionization fraction outside of the completely 
ionized bubbles. 

III. RESULTS 

In order to test the impact of X-ray contribution to reionization on the kSZ effect we 
produce several different simulations using the model above. The model parameters which 
can be adjusted are N ion , C, M min , M ion , X frsiiC , and A. We produce several simulations with 
different parametrizations, but which have an optical depth to the surface of last scattering, 
r, consistent with WMAP measurements (r = 0.088 ± 0.015) [3]. For simplicity, we ignore 
the residual ionization from cosmological recombination (prior to reionization). From these 
models we can evaluate the effect of a varying level of X-ray contribution to reionization. The 
values of the various parameters for each model are listed in Table HI Due to the constraint 
on r, our models have fairly similar global reionization histories. These are shown in Fig. [TJ 

The main difference with X-rays is that their relatively long MFP would cause some 
smooth ionization fraction throughout the IGM in regions which are not inside bubbles. 
This is illustrated clearly in Figure [2j We plot the ionization fraction through a slice of our 
simulation at z = 9 for the "fiducial" and "50% X-ray" models. In both cases we have used 
the same underlying density field and the global mass weighted ionization fraction is roughly 
the same. There is a striking difference in the character of reionization between these two 
models. The UV-only case is dominated by completely ionized bubbles with surrounding 
neutral regions. The model with X-rays has smaller bubbles surrounded by a partially 
ionized IGM. 

In order to quantify the change in the kSZ signal due to X-rays, we compute the angular 
power spectrum of CMB temperature fluctuations, C[. We calculate the kSZ temperature 
fluctuations from the EoR [z = 6 — 20) by integrating Eq. ([1]) through our simulations. The 
resulting 2D temperature field is then Fourier transformed to compute the power spectrum. 
We include the kSZ signal from lower redshift by calculating the power spectrum analytically 
and adding it to the EoR result. This is necessary because at low redshifts non-linear effects 
become significant and they are not accurately captured with the density fields we use 
from 21cmFAST. Wc calculate the analytic power spectrum using the method described 
in [2[ (see their Eq. (4)). To account for the primordial temperature fluctuations, which 
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z 



FIG. 1: The mass weighted global ionization fraction as a function of redshift for the different 
reionization models. 

dominate on large angular scales, we add the lensed C/'s calculated from CAMB [loj]. We 
make the assumption that the thermal Sunyaev-Zel'dovich effect, with its distinctive spectral 
signature, will be removed completely. 

All of the models tested were designed to have nearly identical optical depth due to 
electron scattering. Of course there will always be some uncertainty in value of this quantify. 
For example, Planck measurements will constrain r to within « 0.005 jUJ. We find that this 
uncertainty should not strongly affect our results. As an example, for the "fiducial" model, 
we vary iV ion within the corresponding uncertainty in r and find that the power spectrum 
never goes as low as the 10% X-ray model for the I values in Figure |3j 

In addition to plotting the theoretical estimates, we also show the error bars of SPT 
in Fig. H] [22| . SPT will only be sensitive to the kSZ signal at I w 3000 due to the high 
amplitude of the primordial temperature fluctuations at lower I and dusty star-forming 
galaxies at higher I. Detailed estimates b y th e SPT team project an ultimate sensitivity of 
A (Ql(l + 0/(2tt)) = l/iK 2 /T^ MB . ACT Q will be a couple times less sensitive (Matthew 
McQuinn, private communication 2011). 

IV. DISCUSSION AND CONCLUSIONS 

We have calculated the effect of an X-ray contribution to reionization on the fluctuations 
in the CMB from the kSZ effect. This was accomplished with semi-numerical simulations 
of reionization. For the case with no X-ray contribution, we tested how changes in the 
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FIG. 2: The ionization fraction through a slice of our "fiducial" (left panel) and "50% X-ray" (right 
panel) model at z = 9. Note that the same underlying density field has been used for both. The 
slice corresponds to a 170 x 170 x 0.66 Mpc 3 volume in the simulation. The color bar gives the 
mass weighted neutral fraction in each pixel. 
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5000 10000 
I 

FIG. 3: The angular power spectrum of temperature fluctuations for the different reionization 
models tested. 
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FIG. 4: The angular power spectrum of temperature fluctuations at I = 3000, where ACT and SPT 
will be sensitive to kSZ signal. The expected accuracy of SPT, A (QZ(1 + l)/(2n)) = l/xK 2 /Tj~, MB , 
is shown as an error bar. The various reionization models are shown with the same lines as Fig. [3l 



physical parameters, including the number of ionizing photons produced, the dumpiness of 
the IGM and the minimum mass of dark matter halos which can host galaxies affect the 
kSZ signal. This was then compared with various levels of X-ray contribution to reionization 
parametrized by the fraction of ionizing photons which are X-rays, Xf rac . 

Our reionization model assumes that a constant number of ionizing photons are produced 
for each baryon which is incorporated into galaxies. We then calculate how many baryons 
have collapsed into dark matter halos at a given redshift using the extended Press-Schechter 
formalism. UV photons have a MFP much smaller than the pixel scale in our simulations 
while X-rays are assumed to uniformly ionize regions with a fixed comoving MFP, for sim- 
plicity. Varying the MFP between A = 20 Mpc and A = 50 Mpc yields very similar results, 
suggesting that our conclusions are robust. 

Due to the comparatively long MFP of X-rays, they can partially ionize the regions in 
between the ionized bubbles created by UV photons. For a fixed number of total ionizing 
photons this results in smaller bubbles surrounded by partially ionized IGM. This affects 
the kSZ signal by reducing the "patchy" component of the fluctuations, associated with the 
non-uniform topology of reionization. In this way it is possible to distinguish models which 
have a very similar global reionization history but different X-ray contributions. 

In particular, we find that an X-ray contribution greater than Vf rac = 10% produces 
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a substantially lower power spectrum than our UV-only models constrained to have equal 
total optical depth to electron scattering. This suggests that it may be possible to constrain 
the X-ray contribution to reionization with future kSZ observations. Unfortunately, the 
sensitivity of SPT and ACT may be insufficient to distinguish the models presented in this 
paper. However, a non-detection of the kSZ signal from the EoR may suggest a significant 
X-ray contribution. We find that our X-ray models push the kSZ power from reionization 
below the anticipated sensitivity of SPT and ACT. Because annihilation or decay of dark 
matter particles could imprint a similar uniform ionization component, such a non-detection 
could also be used to constrain the properties of dark matter. 
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